Directional interacting whispering gallery modes in coupled dielectric microdisks 
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We study the optical interaction in a coupled dielectric microdisks by investigating the splitting 
of resonance positions of interacting whispering gallery modes (WGMs) and their pattern change, 
depending on the distance between the microdisks. It is shown that the interaction between the 
WGMs with odd parity about y-axis becomes appreciable at a distance less than a wavelength and 
causes directional emissions of the resulting interacting WGMs. The directionality of the interacting 
WGMs can be understood in terms of an effective boundary deformation in ray dynamical analysis. 
We also discuss about the oscillation of the splitting when the distance is greater than a wavelength. 

PACS numbers: 42.55.Sa, 42.65.Sf 



I. INTRODUCTION 

In spherical and cylindrical dielectric cavities high- 
Q modes are the whispering gallery modes (WGMs) in 
which light rays circulate along the curved inner bound- 
ary of the cavities, reflecting from the boundary with 
an incident angle always greater than the critical angle 
for total internal reflection, thus remaining trapped in- 
side the cavities 0, 01 ■ There are only minute isotropic 
emissions of light caused by evanescent leakage. For the 
applications to optical communication and optoelectric 
circuit, this isotropic emission is not desirable, rather di- 
rectional emission is much more useful and effective [2.J. 

As a simple system for directional emissions, slightly 
deformed microcavities have been proposed, and direc- 
tional emissions, tangential from the boundary points 
with the highest curvature, are achieved. In the ray dy- 
namical viewpoint, as being slightly deformed, some in- 
variant tori are destroyed in the Poincare surface of sec- 
tion (PSOS), but the Kolmogorov-Arnold-Moser (KAM) 
tori still confine the rays supporting the WGMs. In this 
case, the tunneling process, through the lowest dynamical 
barrier, to the critical line for the total internal reflection 
can explain the tangential emissions 2, 3]. 

When the cavity boundary is highly deformed, the 
PSOS shows a global chaotic sea with very small 
integrable regions (islands). In this strong chaotic 
case, the directional emissions can be found in scarred 
resonances!^ and quasiscarred resonances!^. Unlike the 
slightly deformed case, the direction of emission in these 
resonances can be deviated from the tangential, and is 
well explained by the unstable manifold structure near 
the critical line for total internal reflection 6] and the 
Fresnel filtering effect 0j- In addition, there are special 
boundary shapes for generating unidirectional emission, 
spiral-shaped 8] and rounded triangle-shaped 9]. We note 
that the efforts for directional emissions are mainly based 
on the deformation of boundary shapes. 

In this paper we show, through a numerical study on 
the interacting WGMs in a coupled identical disks, that 
a mode-mode interaction can generate directional emis- 



sions. The interaction between two WGMs is parameter- 
ized by the distance between two disks, and it turns out 
that the strength of the interaction between WGMs with 
odd parity about y-axis becomes appreciable at a dis- 
tance less than a wavelength, which are evident from the 
results on variation of resonance positions and patterns. 
In order to explain the resulting directional emission, we 
assume that the circular boundary shapes would be effec- 
tively deformed due to the mode-mode interaction. With 
this assumption the ray dynamical analysis gives a good 
explanation for the degradation of Q-factor and the en- 
hancement of directional emissions. In addition, when 
the WGMs are weakly coupled, the resonance positions 
show an oscillating behavior depending on the distance. 

The paper is organized as follows. In Sec. II we il- 
lustrate our system, i.e., a coupled dielectric microdisks. 
The numerical results for the strongly interacting WGMs 
and a ray dynamical model with an effective deforma- 
tion are presented in Sec. III. The behavior of weakly 
interacting WGMs is discussed in Sec. IV. Finally, we 
summarize results in Sec. V. 



II. COUPLED DIELECTRIC MICRODISKS 

As a simple system for the study on interacting WGMs, 
we take coupled dielectric microdisks. In this system, the 
strength of the interaction can be controlled by adjusting 
the distance between two disks. It is well known that the 
WGMs of a single dielectric microdisk are two-fold de- 
generate due to the circular symmetry, and are classified 
by the angular momentum mode index m and the radial 
mode index Therefore, in the coupled microdisks, 
the WGMs would show an approximate four-fold degen- 
eracy when the distance d between the disks is very large, 
equivalently, the interaction between WGMs is negligi- 
ble. As the distance is getting smaller, the WGM in one 
microdisk starts to know the existence of the WGM in 
the other microdisk, then the system is no longer cir- 
cular symmetric, and the four-fold degenerate resonance 
positions start to split each other, and the degree of the 
splitting measures the strength of the interaction between 
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FIG. 1: The four symmetry classes of the coupled identical 
disks. The former letter indicates the parity on x and the 
latter does the parity on y coordinates. Even (Odd) symmetry 
is marked by dashed (solid) lines. 



the WGMs. 

Figure 1 shows the four symmetry classes of coupled di- 
electric microdisks. The system has two symmetry lines 
and the former letter E{0) is even(odd) if the wave func- 
tion is even(odd) with respect to x — > —x and the latter 
refers to y — > —y. 

In practical calculation of the interacting WGMs in the 
coupled dielectric microdisks, we use the boundary ele- 
ment method (BEM) which is effective when the bound- 
ary is strongly deformed from a circular shape and the 
cavities are coupled jlQ]. In this paper, we focus on the 
TM polarization where both the wavefunction and its 
normal derivative are continuous across the boundary. 
The radius of disks is i?, the distance between two dielec- 
tric disks is d, and riin and Uout are the refractive indices 
inside and outside disks, respectively. We set Uout — 1 
throughout the paper. 



III. STRONGLY INTERACTING WGMS 

In this section we present numerical results on the vari- 
ation of resonance positions and patterns of the strongly 
interacting WGMs, i.e., the case of the short distance, 
d < X. As mentioned before, we expect that the four- 
fold degenerate WGMs would start to split each other 
due to the interaction between WGMs as the distance d 
decreases. As a result of the interaction the directional 
emissions appear in the interacting WGMs with odd par- 
ity about y-axis. We explain this directionality by assum- 
ing an effective deformation of boundary in ray dynam- 
ical analysis. In the practical BEM calculation, we take 
12 elements per a wavelength inside (Ai„ = 27r/7T,i„fc). 



A. Variation of resonance positions and patterns 

Numerical calculation is performed for the WGM of 
mode index {m,l) = (77,1) with ni„ — 1.4. Exact res- 
onance positions of the WGM in an isolated microdisk 
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FIG. 2: (color online) The splitting of the degenerate 
WGM(^Y7,i) due to the mode-mode interaction in the coupled 
disks system. Re(fei?) (a) and lm{kR) (b) are plotted de- 
pending on the distance d between two disks. Note that kR 
converges into the resonance position of the single disk case, 
59.7155— iO. 0039 as d increases. Black circle, red square, green 
diamond, and blue triangle correspond to EE, EO, OE, and 
OO modes, respectively. 



can be obtained from the matching conditions between 
the Bessel function and the Hankel function of the first 
kind which are inner and outer radial solutions of the 
Helmholtz equation, respectively. The exact resonance 
position of T4^GM(77^i) is kR = 59.7136 - i2.5687 x 10"^ 
where k is the vacuum wavenumber. The very small 
value of |Im(fci?)| means high Q factors from the rela- 
tion Q = — Re[fci?]/2Im[fci?]. Unfortunately, it is very 
difficult to get the exact imaginary values from the 
BEM, for example, when we take 12 elements per Ai„, 
the BEM calculation for the isolated microdisk gives 
kR = 59.7155 - i3.9 x lO'^ for WCMf^r^^iy In spite of 
the restriction in the precision of resonance positions, we 
rely on the BEM calculation in analyzing the interacting 
WGMs based on the following reasons. First, the change 
of d{kR) ~ 10^'' does not give any visible variation in 
resonance patterns. Second, the BEM calculation would 
give a correct result when the variation of resonance po- 
sitions exceeds the precision limit, S{kR) ~ 10^"^. 

Figure 2 shows the variation of resonance positions of 
the interacting WGMs. The four fold degenerate state 
starts to split into two groups at ~ 0.8A in Re[A;i?] and 
at d ~ 0.4A in Im[fci?]. Here we can sec that the reso- 
nances with odd parity about axis(Oii^, OO modes) 
show larger variations, indicating that the WGMs in 
the resonances are strongly coupled. We note that the 
Ke[kR] values of OE, OO modes increase with decreas- 
ing d. The increment of Re[/ci?] implies a reduction of ef- 
fective boundary perimeter, and the interaction between 
WGMs in the OE, OO modes is, thus, repulsive. From 
the same argument, we conclude that the interaction in 
the EE, EO modes is weakly attractive. This result will 
be used in determining the effective deformation in the 
ray model in the next subsection. As shown in Fig. 2 
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FIG. 3: (color online) Resonance patterns of the interacting 
WGM(77_i)S when Uin = 1.4 and d = 0.02A. (a) EE mode, 
(b) OE mode. Red-white-blue-dark blue colors indicate high 
to low intensity on a logarithmic scale. 



(b) the Ini[fci?] values of OE, OO modes decrease with 
decreasing d, implying the degradation of Q-factor due to 
the repulsive interaction. Therefore, we can expect that 
the emission of the OE, OO modes would be stronger 
than those of the EE, EO modes. From the viewpoint of 
the effective boundary deformation due to the repulsive 
interaction, we can understand the difference of the onset 
points of the splitting in Re[fci?] and Im[fci?] . In a slightly 
deformed cavity, the rays supporting the WGM are com- 
pletely confined by the KAM tori, resulting no drastic 
reduction of Q-factor. As the cavity is more deformed, 
the KAM tori would be broken and the rays can diffu- 
sively escape along the unstable manifolds, and then the 
Q-factor, equivalently Im[A:i?], would decrease rapidly. 

From the splitting behavior of the resonance positions, 
we can expect that the resonance patterns of EE and 
EO modes would be different from those of OE and OO 
modes. In Fig. 3 the resonance patterns of the interact- 
ing WGM(Tj i)S are shown when the distance d is 0.02A. 
The resonance position of EE mode shown in Fig. 3 (a) 
is kR — 59.6829 — ^0.0040 and, as expected from the small 
absolute value of Im[fci?], the very small evanescent leak- 
age is shown [Tlj|. However, as shown in Fig. 3 (b), the 
resonance pattern of OE mode {kR = 59.7998 - i0.0215) 
shows clear directional emissions, reflecting the strong 
repulsive mode-mode interaction. Its far field emission 
pattern is plotted in Fig. 5 (a) where the four strong 
directional emissions are clearly seen. We find that the 
directionality of the emission pattern is insensitive to the 
distance d, although the strength of emissions decreases 
with increasing d. 

In OE and OO modes, the strong directional emis- 
sion by the repulsive mode-mode interaction is a generic 




FIG. 4: (color online) Resonance patterns of the interacting 
VyGM(29,i)S when run = 2.0 and d = 0.005A. (a) EE mode, 
(b) OE mode. Red-white-blue-dark blue colors indicate high 
to low intensity on a logarithmic scale. 



feature, but the emission directions are closely related 
to the reflective index n™. As an example, the reso- 
nance patterns of EE and OE modes in the interacting 
WGM(^29,i)S, when = 2.0 and d = 0.005A, are shown 
in Fig. 4. As expected, four strong directional emissions 
are shown only in the OE resonance pattern in Fig. 4 
(b) , and the two beams emitted from one disk are almost 
parallel to x-axis. Note that the direction of emission is 
quite different from that of the interacting WGM(yT i-) 
in Fig. 3 (b). The corresponding far field emission pat- 
tern is shown in Fig. 5 (b) where we confirm the two 
directional emissions along x-axis and clear interference 
pattern of the parallel beams. 

In the next subsection we will introduce a ray dynami- 
cal model to explain the Uin dependence of the direction- 
ality of emissions in the interacting WGMs with OE and 
OO parities. 

B. Ray dynamical model : Effective deformation 

In a circular disk, the ray dynamics is simple. The 
rays with incident angles greater than the critical an- 
gle, 9c = arcsin(l/ni„), are perfectly confined in the disk 
by the total internal refiection, while the other rays es- 
cape isotropically due to its rotational symmetry. With 
a simple combination of this trivial ray dynamics, it is 
impossible to explain the directionality of the interacting 
WGMs shown in the previous subsection. 

We recall that the interaction between WGMs are re- 
pulsive in OE and OO modes. In fact, this originates 
from the constraint that the field value at a: = or on 
y-axis should be zero due to the odd parity about y-axis. 
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FIG. 5: The far field emission patterns of the interacting 
WGMs. (a) The interacting WGM(r7,i)S (n™ = 1.4) shown 
in Fig. 3 (b). (b) The interacting WGM^29,i)S {riin = 2.0) 
shown in Fig. 4 (b). 



As a result, the intensity spots confronting each other 
near a; = would shift repulsively, and the structure of 
the whole intensity spots in a WGM would be slightly 
deformed from circle. In order to incorporate this ef- 
fect of the repulsive interaction into the ray dynamics, 
we consider a slightly deformed circular boundary which 
can support the slightly deformed WGM similar to one 
of the interacting WGMs. 

As a simple ray dynamical model containing the effect 
of the repulsive interaction, we consider a circular disk 
with a cut as shown in Fig. 6 (a). The deformation pa- 
rameter is e which is the reduced length of the radius by 
the cut. Figure 6 (b) shows its PSOS representing the 
trajectory of a ray starting from one point in phase space 
(s,p), where s is the boundary coordinate andp = sin0, 9 
being the incident angle, without the consideration of the 
refractive escape. Since the circle with a cut is a discon- 
tinuous deformation(non-KAM system), this model can- 
not describe weakly deformed case where the rays sup- 
porting WGM are still confined in KAM tori. So, our 
model is more suitable to a moderately deformed case 
where the rays supporting WGM can diffusively escape. 
The ray can change their incident angle only through the 
bounce on the cut, and eventually the ray trajectory fill 
up the whole phase space, even though there are so many 
marginally stable lines. The broken lines at p ~ 0.50 and 
p ~ 0.71 represent the families of the marginally stable 
triangular and rectangular periodic orbits, and the gaps 
of lines correspond to the cut of the boundary. 

The dielectric microcavities are open systems where 
rays can refractively escape from the microcavities, 
and the escape rate are determined by the Fresncl 
equations 12]. In order to understand the emission di- 
rection of rays, we obtain the survival probability distri- 
butions in both cases of ni„ = 1.4 and riin = 2.0 which 
are shown in Fig. 7 (a) and (b), respectively. The gray 
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FIG. 6: (a) The boundary shape used in the ray dynamical 
analysis; a circle with a cut. (b) PSOS for the billiard shown 
in (a) when R = 1 and e = 0.01. 



points in Fig. 7 correspond to the rays with normal- 
ized probability greater than 0.1 in the time range of 
50 < t < 53 with a time scale as the length of ray tra- 
jectory, and the rays start from a uniform ensemble of 
1000 X 1000 initial positions in the phase space. The pat- 
tern of the survival probability distribution reveals the 
openness structure on the unstable manifold background 
near the critical lines, Pc — l/f^m, for the total inter- 
nal reflection. The directionality of emissions and the 
emitting part of boundary are explained by the pattern 
below the critical line. The color plots below the critical 
lines in Fig. 7 (a) and (b) illustrate how the long-lived 
rays supporting the WGMs can escape. We find that the 
long-lived rays refractively escape through the unstable 
manifold structure arising in the survival probability dis- 
tribution. As mentioned before, the ray far above the 
critical line changes its angular momentum only when 
bouncing from the cut and, depending on the bouncing 
position on the cut the angular momentum can increases 
or decreases, in other words, the angular momentum dif- 
fuses to other values. From repetition of this diffusion 
process, the ray can reach below the critical line and 
eventually escape the microcavity. The ensemble of the 
above escape process makes the color plots. The dark- 
ness of the color plots represents the population of the 
ensemble. 
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FIG. 7: (color online) The survival probability distribution. 
The colored spikes show the parts through which long-lived 
rays escape, (a) riin = 1-4 case, (b) ni„ = 2.0 case. 



The colored spikes below the critical line(pc — 0.714) in 
Fig. 7(a) correspond to the rays following the diamond- 
typed period orbit inside the microcavity and indicate 
that the ray emission from about s = 37r/2 is very strong 
for the counterclockwise circulating rays(p > 0). These 
are consistent with the resonance pattern shown in Fig. 
3 (b) where the faint diamond structure is seen inside the 
microdisks and the strong emission comes out from about 
s — 37r/2 of the right microdisk. If we consider another 
effective deformed microcavity to simulate the interacting 
WGMs, we can get the resulting emission pattern shown 
in Fig. 8 (a) where we neglect the emission from the cut. 
This ray dynamical result is very similar to the far field 
emission pattern of Fig. 5 (a) except the interference 
oscillation in peaks. The same discussion is valid for the 
case of — 2 case shown in Fig. 7 (b). The resulting 
emission pattern is given in Fig. 8 (b), and this explains 
well the far field pattern of the interacting WGMs shown 
in Fig. 5 (b). 

Although the deformed disk in Fig. 6 (a) is a non- 
KAM system, this explains well the insensibility of the 
directionality of emissions to the distance d. The longer 
distance d corresponds to the smaller e value. In this case, 
although the average escape rate and the emitting part 
of the boundary would decrease, it is clear that the emis- 
sion directions and the emitting positions on the bound- 
ary are essentially invariant. However, the mode-mode 
interaction would create a continuous deformation, i.e., 
the system would be a KAM system. In a KAM model, 
even if the rays supporting the WGM are confined in 



KAM-tori, it is still possible for the rays to reach chaotic 
sea through the dynamical tunneling, and then diffuse 
along the unstable manifolds to the critical linep^. The 
unstable manifold structure near the critical line in the 
KAM model would be similar to that of the non-KAM 
model if the global boundary shapes of both models are 
almost identical. Therefore, both models would give the 
essentially same emission directionality. 



IV. WEAKLY INTERACTING WGMS 

When the distance d is larger than A, the interaction 
between WGMs becomes very small. The strength of 
the small interaction can be measured by the deviation 
of resonance position from that of isolated corresponding 
WGM as done for the strongly interacting WGMs in the 
previous section. Figure 9 (a) shows the variation of the 
resonance position in the range of 2A < d < 6A for the in- 
teracting TyGAf(29,i)S with riin = 2. It is shown that the 
resonance positions of both iJE' (black circle) and OE{red 
square) modes oscillate with the period Ad ~ A/2. The 
oscillatory behavior also appears for EO and OO modes. 
The interference effect of emitted waves from the WGMs 
seems to be responsible for the oscillatory behavior. As 
pointed out in the previous section, the BEM calculation, 
however, has the precision limit of 5{kR) ~ 10~^. We 
note that the amplitude of the oscillations is almost same 
order with the precision limit. So, we have to be careful 
to accept the oscillatory behavior as a real phenomenon, 
because it can be a numerical artifact. To check this, 
we performed the same calculation for the interacting 
T4^GM(7 2)S {riin — 2.0) which are relatively low-Q res- 
onance modes. Since the emission of the WGM(y 2) is 
stronger than the WGM(29,i) case, we expect the varia- 
tion of the resonance position would be larger than the 
precision of BEM. 
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FIG. 9: (color online) The oscillatory behavior of the res- 
onance position of interacting WGMs when d > A. (a) 
WGM(^29,i) case, (b) W^GA/jy 2) case. Black circle and red 
square denote EE and OE modes, respectively. 



The exact resonance position of the WGM^i 2) in an 
isolated circular disk is kR = 6.5806 — iO.1117 and in the 
BEM kR ^ 6.5808 - i0.1127 giving the same precision 
Hniit of 5{kR) ~ 10"^. For the interacting WGM(^j_2)S: 
the variations of resonance positions for EE and OE 
modes are shovifn in Fig. 9 (b). When d > X, both 
real and imaginary parts of kR oscillate with a period 
Ad ~ A and the oscillation amplitude is much greater 
than the precision limit of the BEM, which supports that 
the oscillatory behavior of the interacting WGMi^29,i) in 
Fig. 9 (a) would be a real phenomenon, not a numerical 
artifact. We note that the oscillation of the resonance po- 
sition of OE mode is out of phase with that of EE mode 
in both Re(A;i?) and Im(fci?). The resonance patterns at 
the local maxima of Im(fci?), corresponding to minimum 
leakages of the system, denoted by arrows in Fig. 9 (b) 
are shown in Fig. 10, and these explain why the oscil- 
lations of resonance positions of OE and EE modes are 
out of phase and have a period Ad ~ A. The internal 
patterns of the resonances are almost invariant, and the 
characteristic difference is the number of intensity spots 
on the horizontal axis between two disks which increases 
one by one. Odd and even number of spots appear in OE 
and EE modes, respectively. Roughly we can understand 



FIG. 10: (color online) Resonance patterns of those corre- 
sponding to the maxima of the oscillation indicated by ar- 
rows in Fig. 9 (b). The resonance positions are (a) kR = 
6.5744 - iiO.0712 {OE mode), (b) kR = 6.5824 - i0.0793 {EE 
mode), and (c) kR = 6.5812 - i0.0877 {OE mode), run = 2.0 
is used and red-white-blue-dark blue colors indicate high to 
low intensity on a logarithmic scale. 



this oscillation behavior as the degree of the accordance 
with the quantization condition of the unstable periodic 
orbit lying on the horizontal axis between two disks even 
though the boundary condition on the ends of the period 
orbit is not trivial. In fact this kind of explanation is 
valid only for relatively low kR case where the interfer- 
ence on the unstable orbit dominates, and other inference 
effects are negligible. In the high kR case many beams 
emitted from the intensity spots inside disks take part 
in the interference process, and the resulting oscillation 
of each symmetry mode would show more complicated 
behavior in its period and amplitude. 



V. SUMMARY 

In the coupled disks system, we have shown that the 
strongly interacting WGMs with odd parity about y-axis 
gives good directional emissions, and the directions of 
the emissions are determined by the refractive index of 
the dielectric disks. This finding has been well explained 
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